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Essential spectrum of local 
multi-trace boundary integral operators 

X. Claeys*’* 


Abstract 

Considering pure transmission scattering problems in piecewise constant media, we 
derive an exact analytic formula for the spectrum of the corresponding local multi-trace 
boundary integral operators in the case where the geometrical configuration does not 
involve any junction point and all wave numbers equal. We deduce from this the essential 
spectrum in the case where wave numbers vary. Numerical evidences of these theoretical 
results are also presented. 


Introduction 

Many applications involve the simulation of wave propagation phenomena in media with 
piecewise constant material characteristics that can be modeled by elliptic partial differential 
equations with piecewise constant coefficients. In such situations, the computational domain is 
naturally partitioned into sub-domains corresponding to the values of material characteristics. 

As regards numerical approaches to be used to tackle wave propagation problems, although 
one may opt for finite elements or similar volume methods, boundary integral equation meth¬ 
ods provide accurate alternatives that are less prone to such undesirable effects as numerical 
dispersion. Admittedly discretization of boundary integral equations leads to dense ill con¬ 
ditioned matrices which raises numerical challenges and requires careful implementation, but 
many progresses achieved in the past decade (efficient preconditioners, fast multipole meth¬ 
ods, adapted quadrature techniques) now place boundary integral equation techniques as a 
serious numerical alternative for high performance computations. 

In the context of parallel computing, it becomes desirable to embed integral equation 
methods into a domain decomposition paradigm. The Boundary Element Tearing and In¬ 
terconnecting method (BETI) was developed in this spirit, more than a decade ago, as an 
integral equation counterpart of the FETI method, see [8, 13, 14, 17, 23]. An alternative 
approach dubbed Multi-Trace formalism, leading to different solvers, was introduced a few 
years ago [18, 19, 10, 11, 2, 3, 4, 5], providing other boundary integral formulations adapted 
to multi-domain geometrical configurations. Multi-trace boundary integral formulations seem 
well adapted to block preconditioners for domain decomposition but still very little is known 
about associated iterative global solvers. To our knowledge, the only contributions in this 
direction are [11, 7]. 

A precise knowledge of the spectral structure of the equation under consideration is most 
of the time a key ingredient for devising efficient domain decomposition strategies. The 
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main purpose of the present contribution is thus to provide results concerning the spectrum 
of local multi-trace operators. In particular we describe the spectrum through an explicit 
formula in the case where the propagation medium admits uniform characteristic material 
parameters, which yields the location of the accumulation points of the spectrum of local 
multi-trace operators in the general case of piecewise varying effective wave numbers. In 
addition, the present contribution provides a proof of the well posedness of the local multi¬ 
trace formulations in the case of non-trivial relaxation parameters, which was mentioned in 
[11, Rem.2] as an open problem. 

The present contribution is organized as follows. In the first section we describe the 
scattering problem and the geometrical configurations under consideration. In Section 2 and 
3 we introduce notations related to traces and integral operators, and recall well established 
results of classical potential theory. In Section 4 we recall the derivation of (relaxed) local 
multi-trace formulations as presented in [11], and prove uniqueness of the solutions to this 
formulation for general values of the relaxation parameter. Section 5 presents detailed calculus 
achieved in the case of two particular elementary geometrical configurations. Section 6 is 
dedicated to the study of the spectrum of the local multi-trace operator, and Section 7 will 
present numerical results confirming the theory. 


1 Setting of the problem 


In this section, we will mainly introduce notations, and write the problem under consideration, 
starting from a precise description of the geometrical configurations we wish to examine. 
First of all, we consider a partition M d := U ” =0 flj where the flj’ s are Lipschitz domain, 
where M. d = M,M 2 or M 3 . We assume that each flj is bounded except fig. Changing the 
numbering of sub-domains if necessary, we may assume without loss of generality that each 
Qj is connected. We shall refer to the boundary of each sub-domain by Tj := dflj, and also 
set r j k := T ,j nTj. = dflj n dflk to refer to interfaces. The union of all interfaces will be 
denoted 


n 

E := U r ? - = U ] 

j =0 0<j<k<n 


j,k 


This set will be referred to as the skeleton of the partition. We make a further strong geometric 
hypothesis, assuming that the sub-domain partition under consideration does not involve any 
junction point, so that each Tj ^ is a closed Lipschitz manifold without boundary, 


dTj k = 0 Vj, k = 0 ... n. 


( 1 ) 



fig 
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Sobolev spaces We need to introduce a few usual notations related to standard Sobolev 
spaces. If to C is any Lipschitz domain, we set H 1 ^) := {u £ L 2 (cu), Vu £ L 2 (cu)} equipped 
with the norm ||u|| 2 1(a;) : = + II Vv IIl 2 (u,)’ and Hl ( A > w ) : = € H 1 ^), Av £ L 2 (w)} 

equipped with II' U IIh 1 (A,oj) := IMIh^w) + H Av lliW In addition, if H(w) refers to any of 
the above mentioned spaces, then Hi oc (aJ) will refer to the space of functions v such that 
ipv £ H(w) for any i\) £ < ^ 3 mp (K d ) := {p £ ^°°(M d ), supp(<^>) bounded}. 

For any Lipschitz open set cj C M d , we shall refer to the space of Dirichlet traces 
H i/ 2 (<9w) := {v| 9tJ , v £ H 1 ^)} equipped with the norm |M| H i/ 2 (&y) := min{||«|| H i( w ), u\ dui = 
v}. The space of Neumann traces H~ 1 // 2 ((9w) will be defined as the dual to H l / 2 {duj) equipped 
with the corresponding canonical dual norm \\p\\ u -i/ 2 {du}) := sup. ugH i / 2 (ai;) y {0} |(p, 1 / 1 |^|| H i/ 2 (acv) 

Transmission problem We will consider a very standard wave scattering problem (so- 
called transmission problem), imposing Hemholtz equation in each sub-domain, as well as 
transmission conditions across interfaces: find u £ Hj (M d ) such that 

—A u — k 2 u = 0 in Qj Vj = 0... n 
u — u; nc is outgoing in fig 

< ( 2 ) 
u kj ~ u \r k = 0 

drijUlvj + d nk u\r k =0 on T j:k = Tj n T k , Vj, k = 0... n 

In the equation above ui nc £ H 1 1 oc (M d ) is a known source term of the problem satisfying 
—Artinc — KgUinc = 0 in M rf . In addition, we assume that kj > 0 for all j = 0... n. Problem 
(2) is known to admit a unique solution, see [20] for example. The outgoing condition in (2) 
refers to Sommerfeld’s radiation condition, see [15]. A function v € H 1 1 oc (A,I2o) will be said 
K-outgoing radiating if lim^oo f &B \d r v — inv\ 2 da = 0 where B. r refers to the ball of radius r 
and center 0 , and d r is the partialderivative with respect to the radial variable r = \x\. 

Trace operators As this problem involves transmission conditions, and since we are inter¬ 
ested in boundary integral formulations of it, we need to introduce suitable trace operators. 
According to [22, Thm. 2.6.8 and Thm 2.7.7], every sub-domain Qj gives rise to continuous 
boundary trace operators 7 ^ : Hj^fh,) —> H 1 / 2 (dVlj) and 7 ^ : H | l oc (A, i\ 3 ) —> H - 1 / 2 (3Jlj) 
(so-called Dirichlet and Neumann traces) uniquely defined by 

li{p) := (f\ 9 Q j and 'yi(p) := • \7p\gfy Vp € ^°°(%)- 

In the definition above nj refers to the vector field normal to dttj pointing toward the exterior 
of Qj. We will also need a notation to refer to an operator gathering both traces in a single 
array 

7 J 0) := (7d(u),Tn(u)) 

We shall also refer to 7 d c > 7 n,c defined in the same manner as 7 d, 7 n with traces taken from 
the exterior of Qj. In addition, we set 7 c(u) := {pin, c { v )i 7n,c( v ))- We will refer to mean values 
and jumps to these trace operators, setting 

W («)}:= \ ( l 3 («) + 7c («)) and W («)] : = l 3 («) ~ 7c («) 
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2 Trace spaces 


We want to recast Problem (2) into variational boundary integral equations set in trace spaces 
adapted to the present multi-sub-domain context. The most simple space we can introduce 
consists in the multi-trace space [3, Sect. 2.1] i.e. the cartesian product of local traces: 

H(E) := H(T 0 ) x ••• x M(T n ) 

i i (3) 

where H(T i ) := H + 5(r i ) x H"^) . 

We endow each HI(Tj) with the norm ||(u, ll] H r(r J ) := (II'^IIh 1 / 2 ^■)!/ 2 (rq) 1 ^ 2 ’ an d equip 
H(E) with the norm naturally associated with the cartesian product 

II u IIh(E) := ( II u oIIh(To) ^ + ll u «lliii(T n ) ) 2 

for u = (uo,..., u„) € H(E)t In the sequel we shall repeatedly refer to the continuous operator 
7 : n j=o H ioc( A ^i) H ( s ) defined by y(u) := ( 70 H, ■ ■ ■ ,7n(«)), where n^H^A, %) 
should be understood as the set of u G L 2 oc (M d ) such that u\nj € H^fA, ilj) for all j. We 
also need a bilinear duality pairing for H(Tj) and H(E); writing (, )p. for the duality pairing 
between H 1 / 2 (Tj) and H" 1 / 2 (r j ), we opt for the skew-symmetric bilinear form 

n 

[u,*>J := [Uj.Ojlr,- where 

3=0 (4) 

(ft - (ft.ft)ft ■ 

Next, as in [3, Sect. 2.2], [5, Sect. 3.1], we introduce the so-called single-trace space that 
consists in collections of traces complying with transmission conditions. This space can be 
defined by 

X(£) := c1os h( e){ 7(«) = (7>))?=o I « G tf°°(R d ) } (5) 

where the symbol cIosh(e) refers to the closure with respect to the norm on H(E). By 
construction, this is a closed subspace of H(S). Note also that a function u G H^fA, fi 0 ) x 
• • • x H] oc (A, Cl n ) satishes the transmission conditions of (2), if and only if 7 (u) = ( 7 J (u ))” =0 G 
X(E). In particular, if u G Hj^A, R d ) then 7 (it) = ( 7 J ’(u ))” =0 G X(T). In the sequel, we 
will use this space to enforce transmission conditions. The single-trace space admits a simple 
weak characterization, see [3, Prop.2.1]. 

Lemma 2.1. 

For any u G H(E) we have u € X(E) ][u, 0 ] = 0 Vo G X(E). 

3 Summary of potential theory 

In this paragraph, we shall remind the reader of well established results concerning the inte¬ 
gral representation of solutions to homogeneous Helmholtz equation in Lipshitz domains. A 

* Functions in Dirichlet trace spaces like H 1//2 (9Pj) will be denoted by u,v,w, whereas we use p,q,r for 
Neumann traces. Small fraktur font symbols u, 0, ra are reserved for Cauchy traces, with an integer subscript 
indicating restriction to a sub-domain boundary. 
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detailed proof of the statements contained in the present paragraph can be found for example 
in [ 22 , Chap.3]. 


Let & K (x) refer to the outgoing Green’s kernel associated to the Helmholtz operator —A — k 2 . 
For example & K (x) = exp(zAc|o?|)/( 47 r|a?|) in R 3 . For each sub-domain flj , any (v,q) € H(Fj) 
and any x GM. d \ Tj, define the potential operator 


G i(v,q)(x):= f q{y) < g K (x-y)+v(y)Yi j (y)-(V < g K )(x-y)da(y) . (6) 

Jr, 

The operator G f maps continuously H(rj) into H 1 (A,Hj) x H 1 (A,M d \ flj), see [22, Thm 
3.1.16]. In particular G J K can be applied to a pair of traces of the form D = 7 J ( v ). This 
potential operator can be used to write a representation formula for solution to homogeneous 
Helmholtz equations, see [22, Thm 3.1.6]. 

Proposition 3.1. 

Let u € H 1 1 oc (Hj) satisfy A u + k 2 u = 0 in flj. Assume in addition that u is kq- outgoing if 
j = 0. We have the representation formula 


G£ .( 7 J (u)){x) 


u(x) for x € flj 
0 for x eM. d \ flj . 


(7) 


Similarly, ifv€ H 1 1 0 C (M d \Hj) satisfies A v+k 2 v = 0 in M d \flj ) and is Kj -outgoing radiating in 
the case j 0, then we have G J Kj ('yi(v))(x) = —v(x) for x € M d \flj, and G 3 Kj { r Y > c (y)){x) = 0 
for x € flj. 

The potential operator ( 6 ) also satisfies a remarkable identity, known as jump formula, de¬ 
scribing the behavior of G 3 Kj (v)(x) as x crosses Tj, namely 

M'Gi» = o Vt)GH (Tj) (8) 


which also writes [y- 7 ] • G J Kj = Id. We refer the reader to [22, Thm.3.3.1] (the jump formulas 
are more commonly written in the form of four equations in the literature). Proposition 3.1 
shows that, if u is solution to a homogeneous Helmholtz equation in flj (and is outgoing if 
j = 0) then y 7 • G J Kj ('y 3 (u)) = 7 J (u). This actually turns out to be a characterization of traces 
of solutions to homogeneous Helmholtz equation. 


Proposition 3.2. 

Define G l f(fj) := {7 7 (it) | u € H 1 l oc (H J ) , A u + n 2 u = 0 in flj, u outgoing if j = 0 }. Then 
7 7 • Gi : H(Fj) —>• M(Tj) is a continuous projector, called Calderon projector interior to flj, 
whose range coincides with S“(%) i.e. for any 0 € H(Fj) 

i j • Gi(x,) = t) « t) g e; :(fij ). 

Similarly, defining C° ut (flj) := {7 i(u) \ u € Hj L oc (Oy), Au+k 2 u = 0 inM. d \flj, u outgoing if j 
0 }, we have 7 7 • g£(d) = 0 if and only if D G G° ut (Qj). 

For a detailed proof, see [22, Prop.3. 6 .2]. We shall repeatedly use this characterization as a 
convenient way to express wave equations in the sub-domains flj. Here is another character¬ 
ization of the space of Cauchy data which was established in [3, Lemma 6.2]. 
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Lemma 3.1. 

Consider any j = 0, ... n, and any k £ C \ {0} such that > 0, > 0. Then 

for any u £ E^r.,-) we have: u £ C“(IL,) [u, o]r = 0 Vd £ C™(%). Similarly we have 

u £ C° ut (%) <=► [u, x>] Tj = 0 Vo £ C° ut (%). 


The results that we have stated above hold for any j = 0 ...n. We also set C Q (£) := 
C“ 0 (Oo) x ••• x C“ n (O n ) for a = in,out. The notations just introduced allow a condensed 
reformulation of the well-posedness of (2), see [3, Prop.6.1] for a detailed discussion and proof. 

Lemma 3.2. 

X(E) © e in (£) = H(E). 


To handle Calderon projectors in a multi-sub-domain context, it is more comfortable to 
introduce global operators, so as to reduce notations. First we introduce the continuous 
operator A : EI(E) —»• H(E) defined by 


[A(u),t»] := E"=.o [ Ai (Uj),t)i)]r J . 

with A j := 2 { 7 J } • G j Kj 


(9) 


For all u = (Uj)” =0 , D = (t)j)” =Q £ H(E). Observe that (Id ± A)/2 are projectors, according 
to Proposition 3.2, with ker(A — Id) = range(A + Id) = C m (X) and ker(A + Id) = range(A — 
Id) = C out (E). With this notation, a direct consequence of Proposition 3.2 is that (A - 7 ) 2 = 
( 2 { 7 j } • Gi ) 2 = Id. This is summarized in the next lemma. 


Lemma 3.3. 

(A ) 2 = Id. 

Remark 3.1. In the case of two sub-domains M d = Oo U Hi and one interface E = To = Ti, 
with k = no = ki, there is a remarkable identity relating A 0 to A 1 . Indeed, in this situation, 
the only difference in the definition of G° and G* comes from no = — ni. In particidar, 
denoting 


we have G°(u) = —G \(Q • u) Vu £ H 1 /, 2 (ro) x H 1 / 2 (ro). Note that {7 1 } = Q ■ {7 0 }, so 
multiplying the previous equality by {7 1 } yields Q ■ A 0 = —A 1 ■ Q. □ 

Lemma 3.3 above shows directly that A is invertible. It also satisfies a generalized Garding 
inequality. The next result is proved for example in [20, §4]. 

Proposition 3.3. 

Define 0 : H(E) ->• H(E) by 0(u) := (Q(u 0 ), Q(ui),..., Q(u n )) for all u = (u 0 ,ui,... ,u n ), 
where Q(u,p ) := ( u,—p ) for any ( u,p ) £ H 1 / 2 (r :( ) x H _ 1 / 2 (Fj) and any j. There exists a 
compact operator K : H(E) —> H(E) and a constant C > 0 such that 

$fte{-[Au, 0(u)] + [Ku, u] } > C ||u||| (E) Vu £ H(E). 
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4 Local multi-trace formulation 


In this section we would like to recall the derivation of the local multi-trace formulation in¬ 
troduced in [11], and provide detailed analysis for it. The formulation considered in [10] is 
a particular case of the formulations introduced in [11], But only an analysis of the formu¬ 
lation of [10] has been provided so far. In addition, we should underline our hypothesis (1) 
that discards any junction point in the sub-domain partition, while [10, 11] also considered 
geometrical configurations with junction points. 

A key ingredient of local multi-trace theory is an operator yielding a characterization of 
transmission conditions of (2). Considering u = (uk,Pk)k=o an d 11 = ( v j^Qj)‘j= o> we define the 
transmission operator II : H(£) —>• H(£) by 

{ Vj = + u k 

on T ])k \/j,k = 0...n, j^k. (10) 

Qj = ~Pk 

Clearly, for any function u € H^^A, M d ) we have 7 (74) := ( 7 ■*(u))'j =0 £ H(X) with 7 (u) = 
11 ( 7 ( 74 )). Conversely, considering any function u £ L 2 oc (M d ) such that 74 ^ £ H 1 1 oc (A, flj) for 
all j = 0 ... n, then 7 (it) € H(S) is well defined, and if 7 (u) = 11 ( 7 ( 74 )) then u € Hj^A, R d ). 
Routine calculus shows the following remarkable identities 

n 2 = id, n(o) = n(o) and [n(u),t»] = pi(o),uj Vu,bgi(s). (ii) 

As is readily checked, the operator II maps continuously H(X) onto H(S) under Assumption 
(1) that each T h y. is a Lipschitz manifold without boundary. Elementary arguments on trace 
spaces show that, for any u £ H(X), we have u £ X(S) u = II(u). Since II 2 = Id, this 

can be simply rewritten in the following manner. 

Lemma 4.1. 

range(II + Id) = X(X). 

Now consider u = ( 7 ^ (tx))"_o the traces of the unique solution u to Problem (2). The homo¬ 
geneous wave equation satisfied by u in each sub-domain can be reformulated by means of 
Calderon projectors (A — Id)(u — Ui nc ) = 0. Choose any relaxation parameter a £ C \ {0} 
and add a (Id — II)u = 0 to this equation, which is consistent since u = 7 ( 74 ) must satisfy the 
transmission conditions of (2). Denoting f := (A — Id)iii nc , what precedes implies 

J u £ H(X) and 

\ [ (A — n Q )u, 0 ] = [f, 0 ] Vt> £ H(X). ( 12 ) 

where II a := (1 — a)Id + all 

Observe that the operator of the formulation above can also be rewritten in the form of a 
convex combination A — II a = (1 — a)(A — Id) + a(A — II). Existence and uniqueness of 
the solution to this formulation has already been established only in the case a = 1, see [10, 
Thrn.9 and Thm.ll], For all other values of a, well-posedness of this formulation was an open 
problem so far, as mentioned in [11, Rem.2]. Below we prove uniqueness of the solution to 
(12) for any a € C \ {0}. 
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Proposition 4.1. 

Ker(A - n Q ) = {0} VaeC\{0}. 


Proof: 

Take u = (u,-)™ =0 satisfying (A — Id)u + a (Id — II)u = 0. Thus we have ro := a (II — Id)u = 
(A —Id)u G range(A —Id)nrange(II —Id). Denote t Vj the component of tu associated to dQj so 
that ro = (roo, • • • , tu n ), and set ipj{x) := G J Kj («).,•)(*). We have (A + Id)tu = (A 2 — Id)u = 0, 
which can be rewritten y- 7 • G«-(t Vj) = ^(ipj) = 0 for all j = 0... n. Since we also have, by 
construction, —A ipj—n 2 ipj = 0 in f lj, we conclude that ipj = 0 in Qj, and t Vj = (tVj) = 

[7 J (Vh')] = — 7c(Vj)- If we can pi'ove that ipj = 0 in M d \ Qj for each j, this will show that 
lu = 0. 


Now observe that, since ro G range(II — Id) we have n(tu) + ro = 0 i.e II(tu) = — tt>. As a 
consequence the functions ipj satisfy an homogeneous problem with ’’anti-transmission con¬ 
ditions” 


—A ipj — Kjipj = 0 in M“ \ Qj Vj = 0 ... n , 


ijipj =0 in 
ipj is outgoing (with respect to Kj) for j / 0 
+ 7d , c (^fc) = 0 and 

[ yL,Mj) - = 0 on Tj nr k Vj, k. 


(13) 


Since II 2 = Id and II(iu) = —ro, we have 2ro = lu — II(tu), and ro + II(ro) = 0. From this and 
(11), we obtain 2[lt), ro] = [ro — II(tu), tu] = —[to + n(tn), tti| = 0. This can be rewritten 


n p 

/ TD,c(V’j)7N,c(V’j)^ } = 0. 


1=0 


(14) 


Take r > 0 sufficiently large to garantee that M. d \ Do C B r where B r C M. d refers to the ball 
centered at 0 with radius r. Since — Aipj — K 2 ipj = 0 in M d \ Qj, applying Green’s formula in 
each B r \ Dj yields 


Id B r 


ipjd r ipjda = \Vipj\ 2 - Kj\ipj\ 2 dx + 


iB r \n 


dfli 


0 = 


> B r \Qo 


\Vip 0 \ 2 - nl\ip 0 \ 2 dx + / 


'an 0 


In these equations, ”<9 r ” refers to the radial derivative. Take the imaginary part of the 
identities above, and sum over j = 0... n, taking account of (14). This leads to 


E 1 

l=o 


dB r 


n n 1 

ipjdripjda | = Am | ^ 7 d = 0 


where 1 refers to the imaginary unit. By construction, the functions ipj are outgoing radiating, 
so that 0 = linv^oo f dB \d r ipj — mjipj\ 2 da = 0. As a consequence, we finally obtain 




This shows in particular that lim^oo f dB \ijjj\ 2 da = 0 for all j = 1... n. As a consequence we 
can apply Rellich’s lemma, see [6, Lemma 3.11], which implies that ipj = 0 in the unbounded 
connected component of each M d \ Qj. 

Let us show that tpj also vanishes in bounded connected components of W l \ Clj. Take 
an arbitrary j, and let 0 be a bounded connected component of M. d \ 11,. We have 80 = 
dflj fl 8Qk = rfor some k = 0 .. .n,k ^ j. Let O' be the unbounded connected component 
of M. d \ f\,. Then we have Qj C O', and 80' = 80 = 8Qj n Since -0^ = 0 in O', we have 
7 D,c(V , ,j)|rj )fc = ~ 7 D,c(^fc) = 0 and = 7 n. cW’fc) = 0, according to the transmission 

conditions of (13). Finally we have —A ifrj — n 2r tpj = 0 in 0 with 7 c(Vb) = 0 011 90. We 
conclude by unique continuation principle (see [15, §.4.3]) that if)j = 0 in 0. We have just 
proved that 

V’j = 0 in \Qj Vj = 0 ... n. 


Finally we have, on the one hand, tt> = (A — Id)u = 0, so u = (A + Id)u/2 G range(A + Id) = 
C in (S), and on the other hand XV = a(II —Id)u = 0, so u = (II+Id)u/2 G range(II+Id) = X(S). 
So we conclude that u € C in (E) fl X(X) = {0} according to Lemma 3.2. □ 


5 Examples 

Before going further into the analysis of the local multi-trace formulation (12), we dedicate 
this section to deriving and studying it in ultra simplified situations where all calculations can 
be conducted quasi-explicitly. Here we will systematically consider the case where all wave 
numbers are equal 

«o = «i = • • • = «n- (15) 

This assumption will allow substantial simplifications. Another purpose of the present section 
is to determine the spectrum of the multi-trace operator in these simplified situations. 


5.1 Two domain configuration 

We start by considering the case where the space is partitioned in two domains only. This 
simple case was already considered in [10, §3.1], but here we are going to formulate additional 
comments. In this case E = To = Ti = Top. We want to represent the operator (1 — a)(A — 
Id) + a(A — n) : H(ro,i) x IHI(r0,1) —> E^Top) x lHI(r0,1 ) in a matrix form. First of all, note 
that the operator n admits the following expression, 



0 Q " ’ d 0 " 

Q 0 di 


with Q : = 



(16) 
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Hence, denoting u = (uo,Ui), we have 

uo 
Ul . 

To determine the spectrum of the (relaxed) multi-trace operator (1 — a)(A — Id) + a(A — n), it 
suffices to determine the spectrum of A—an. If we compute the square of this operator, taking 
account of (15), we obtain (A —an) 2 = A 2 + a 2 n 2 — a(nA +An) = (l + a 2 )Id —a(nA + An). 
In this case, a direct calculus shows that QA° = — A X Q. As a consequence, an explicit calculus 
yields 


(A — «n)u = 


A 0 -aQ 
—aQ A 1 


-1 

> 

0 

0 


' 0 Q ' 


O 

> 

O 

<0 

_l 

0 A 1 


Q 0 


A X Q 0 j 


0 

1 

.O 

> 

_ 1 


' 0 Q ' 


'A 0 O' 

QA° 0 


Q 0 


0 A 1 


From this we conclude that nA + An = 0, which finally yields (A — an) 2 = (1 + a 2 ) Id. This 
expression, together with the spectral mapping theorem ( see [21, Thm. 10.28] for example), 
provides an explicit characterization of the spectrum in the case where all wave-numbers are 
equaF 

©(A — n a ) C { — 1 + a + \A + a 2 , —1 + a — \/l + a 2 }. 

5.2 Three domain configuration 

Now we consider a partition in three domains M. d = Ho U fR U fR. This situation is pictured 
below. 



As the definition of the transmission operator is given by a formula written on each inter¬ 
face, let us decompose traces on each sub-domain according to interfaces. Considering the 
decomposition To = Top U To, 2 j any trace n G H(To) induces an element R?(t)) € IHI(ro,i) de¬ 
fined by Ri(o) = t)|r 01 - We may define R^o) G lHI(r 0 , 2 ) similarly. This establishes a natural 
isomorphism (R^R®) : H(To) —> IHI(r 0 , 1 ) x m(To, 2 )- The adjoint of those maps are extension 
operators i.e. (R°)*(t>) = t> • 1 r 0 :I G H(To) for any D G H(roj). With these maps, the operator 

5 In the remaining of this article, the square root of complex numbers shall be defined by \/pexp(i6 ) = 
x /pexp(i0/2) for 9 £ (—7t,7t]. 
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A 0 induces a 2 x 2 matrix denoted [A 0 ] with integral operator entries 

a ° A° 

A l,l a 1,2 

A 0 A 0 

a 2,1 a 2,2 


[A°] = 


with A£ fc := R° • A 0 • (RjJ) 


i0 \* 


Plugging this decomposition into the definition of the operator A yields a 4 x 4 matrix of 
integral operators acting on tuples (uo,i, Uo, 2 ,Ui,U 2 ) € H(ro,i) x IHI(r 0 , 2 ) x IHI(r 1 ) x IHI(r 2 ), and 
given by the following formula 


(A — o4I)u = 


A?,! A? i2 -aQ 0 - 


" u 0 ,i ' 

A 2 ,l A 2,2 0 —aQ 


Uo ,2 

-aQ 0 A 1 0 


Ui 

0 -aQ 0 A 2 


. U 2 . 


(17) 


As in the previous paragraph, let us compute the spectrum of the local multi-trace operator 
(1 — a)(A — Id) + a(A — II). Here again, it suffices to determine the spectrum of (A — an). 
Once again we have (A — an ) 2 = A 2 + a 2 n 2 — a(nA + An) = (1 + a 2 ) Id — a(nA + An). 
Besides, taking account of (15), a direct and thorough calculus shows that Q • Ay • • Q = —AT 
So if we compute explicitly the expression of nA + An taking account of this identity, we 
obtain 


nA +An 


0 0 0 A? 2 Q ■ 

0 0 A^Q 0 

0 qa ? 2 0 0 

QA° a 0 0 0 


This time we have nA + An 7 ^ 0. Let us compute the square of this operator. Since Q 2 = Id, 
we obtain 


(nA + An) 2 


A 1,2 A 2,1 0 

0 A 2,1 A 1,2 
0 0 

0 0 


0 0 

0 0 

Q A 1/2 A 2,1 Q 0 

0 QA°iA° 2 Q _ 


Now let us have a closer look at each of the operators A ) 1 2 A 2 x and A 2 1 A) ) 2 . Observe that 
"AJ* = 2 Q • • G k ■ Q for j / k. According to the second equality in (7), we have 

7 1 • G 2 • y 2 • G 1 = 0. As a consequence, we obtain 

A?, 2 A^ = 4 (Q • 7 1 • G 2 • Q) • (Q • 7 2 • G 1 • Q) 

= 4 Q • ( 7 1 • G 2 • 7 2 • G 1 ) • Q = 0 


We show in a similar manner that A 21 A5 2 . To conclude we have (nA + An ) 2 = 0. Such 
a nilpotent operator has a non-empty spectrum (see [21, Thm.10.13]) that is reduced to 
{0} according to the spectral mapping theorem [21, Thm.10.28], which implies of course that 
6 (nA + An) = {0}. Finally we obtain the following spectrum, like in the previous paragraph, 

6 (A — n a ) C { — 1 + a + y/l + a 2 , —1 + a — yjl + a 2 }. 
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6 Spectrum of the operator in a general configuration 

For both examples of the previous two paragraphs, the spectrum of the local multi-trace 
operator only consisted in the two eigenvalues — 1 + a ± \/l + a 2 in the case where all wave 
numbers equal. Besides, during the calculations above, the geometry of the interfaces never 
came into play. In the present section we will show that these are actually general results 
that hold for any number of sub-domain arranged arbitrarily, provided the geometry does not 
involve any junction point. 


To investigate this question in the general case, we need to introduce further notations. 
Recall that each boundary can be decomposed in the following manner F ? = U^jTj^. We 
will decompose traces accordingly. For a given pair j. k with j / fc we define 

R { : H(r i ) H(r i)fc ) 

R £(f) := fir,.. 

To reformulate the above definition fully explicitly, for n = (v,q) £ IH^r^) = H 1 / 2 ^^*,) x 
H- 1 / 2 (r jifc ) and u = (u,p) £ HfFj) = H 1 // 2 (r,,) x H -1 / 2 ^,,), we have 0 = R^,(u) if and only if 
u|r, k = v and p |r, k — Q- The adjoint of these restriction operators are given by 


(Rfc)* : H(F . jlk ) H(r,-) 

(Rfc)*(f) :=f'lr,, fc 

Decomposing traces with the embedding/restriction operators that we have just defined, each 
A 3 : H(F ; ) —> H(r i ) induces a matrix of integral operators denoted [A J ] with maximal size 
(n — 1) x (n — 1) given by 


[A J ] 





A j 


:= R i • V ■ (Rj 


( 18 ) 


In the notation above, it should be understood that the rows and the columns associated to 
indices k such that F ? y 0 = dQj (~l dQk = 0 must be omitted. The row/column associated to 
k = j is to be omitted as well. For each sub-domain Qj, we will need to consider the set of 
indices 

J?j := { k | k ^ j, dflj n / 0}. 

Hence the matrix in (18) is square with card(j^) rows. If we plug the definition of the 
restriction/embedding operators R{ into the definition (9) of the operator A J and the potential 
operator (6), we obtain an explicit formula for Aj, in the case where fc / m namely, 


a u;)-( 




y)[ (n(») x Vu(x)) ■ (n(y) x Vv(y)) 

-K 2 n(x) ■ n {y)u(x)v(y) - p(x)q{y) ] dcr(y)da(x) 


( 19 ) 



(VGrtX* 


y) ■ [ q{y)u(x)n{x) - p(x)v(y)n(y) ] da{y)da(x) 
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for any (u,p) G H(rj >m ) and (v, q) G E^T^fc). Note that the expression above does not hold 
for k = rn. Expression (19) clearly shows that the operators A), m are compact for k rn 
since it only involves smooth kernels. The following lemma yields several remarkable identities 
satisfied by the elements of the decomposition (18). 

Lemma 6.1. 

For any j = 0 ... n, and any k, l , m G J?j with k / m and k / l we have 

*) ( A {,k ) 2 = ld > 

a) Q • A j k k = -Aj j • Q if Kj = K k 

™) A {k ' A l,m = 

Proof: 

Pick an arbitrary j = 0 ... n that will be fixed until the end of the proof. Take a arbitrary 
k, m G J?j with k / m. Let 0 C be the maximal open set satisfying 80 = T^m and 
flj fl 0 = 0. Take an arbitrary 0 G M(Tj )TO ), and denote ti := (R^)*t). Observe that G£(6) G 
H^ oc (M d \ 0). Since T hk n T j jTn = 0, this implies in particular that G«(6) does not admit any 
jump across F^. As a consequence, we have 

A l,J») = 2Ri • l j • G 7 (6) = 2R( • 7 7 • G 7 (6) 

Now observe that, if to G C° ut (f2j), then (R^)*Rj£(tt>) G C° ut (Dj) for any k G J**). Taking m = 
7c -G 7 (Ii), we see that th = (R|,)*A^. m (t>) G G° ut (Qj ). According to the integral representation 
Theorem 3.1, this implies G J K (tt))(x) = 0 for x G f lj. In particular we have 

2Rjyoi(m) = Ajj • Aj m (o) = 0. 

This establishes in). Now we know that (A-') 2 = Id. This implies that, for any k G J?j, we 
have Emef A k m ' A m k = Id- But according to in), all the terms of this sum vanish, except 
for k = m which establishes i). 

To prove ii), observe that Q • R). • {y- 7 } = • {7 fc }. On the other hand, the explicit expression 

of potential operators given by (6) shows that G J K ■ (R),)* = —G(f • (R^)* • Q in the case where 
Kj = K k - Combining these two identities we obtain Q ' A), = 2 Q • R). • {y 7 } • Gf ■ (R).)* = 

-2 Rj • { 7 fc } • G{i • (RJ)* • Q = -Ajj • Q. 

□ 

Next we need to introduce an operator involving only the diagonal blocks of the matrix 
representing A in the decomposition (18), without any term coupling different interfaces. 
Define D : H(E) —> H(E) by the explicit formula 

n 

[D(u), t)J := J2 I Hk Jr,- k Vu, o G H(E) 

j=o keSj (20) 

where u J>k := H) iu ( ). V jtk := R^.(t)j) 


13 



Lemma 6.2. 

We have (D) 2 = Id. Moreover, if all wave numbers equal i.e. «o = «i = • • ■ = n n , we have 

n • d + d • n = o. 

Proof: 

According to Lemma 6.1, we already know that (A^ k ) 2 = Id. Pick an arbitrary pair 

of traces u, 0 G H(£), and set Uj.k := R].(uj) and v J: k := R{,(tr/) for all j, and all k € J?j. 
Applying Formula (20) twice yields 

n 

P», D | := Y, E [(<*) V. 

1=0 k£J?j 

n n 

= E E = Ei u t Mp = i u ’ °i- 

1=0 1=0 

which establishes that D 2 = Id. To establish the second statement, let us first point out that 
the definition (10) of the operator II can be rewritten [II(u), oj = J2j =o Thk&jtj [Q( u fc,l)> Oj.fclr, k 
Combining this with the definition of the operator D, and using Property in) of Lemma 6.1, 
we obtain 


|[n-D(u),o] :=EElQ A b%-* 
j=o keSj 

n 

■= E E - [ A fc,fcQ( u fcj)>» , j,fclr J -,fc = -P • n(u),t»] 

l=o keSj 

□ 

Lemma 6.3. 

The operator T := A — D is compact and satisfies T 2 = 0. 

Proof: 

The operator T only involves terms A). with k fi m. Since P^ n T j jjn = 0, these 
operators defined by (19) only involve smooth kernels. So each A) k ^ m is compact, and 
T is compact itself. Let us compute explicitly the expression of T 2 . Pick arbitrary u,C6 H(S), 
and set u J: p := Rj^Uj) and Ojj- := R^Oj) for all j , and all k G We have 

n 

j= 0 
n 

E E D l,fcl r j,fc 

j=0 

To conclude it remains to apply property ri) of Lemma 6.1. Since, for each term A{ ,Aj we 
have k l and m l , the whole sum vanishes, and we have [T 2 (u), o] = 0. □ 


[T(u),t>] : = 

[T 2 (u),ofl : = 
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Note that it is a direct consequence of the above lemma and Proposition 3.3 that D also 
satisfies a Garding inequality. In addition, since II 2 = Id, a recurrence argument combined 
with the previous lemma readily leads to the following corollary. Note that, remarkably, this 
result holds without any assumption on the wave numbers. 

Corollary 6.1. 

We have (FIT + TII) fc = (IIT) fc + (Tn) fc , \/k > 0. 

Now let us formulate a few elementary and useful remarks concerning the geometrical ar¬ 
rangement of the interfaces. Let T = {0,1,..., n }, and say that two indices j. k are adjacent 
if dQj n 0. This adjacency relation endow T with a graph structure. We order the 

elements of T, writing j -< k if Clj is included in a bounded connected component of \ I\,. 
This induces an tree struture on T. In particular T is a tree and does not admit any (simple 
cycle), see for example [1, Chap. 16]. And it does not admit chain with a length larger than 
n. The picture below provides an example of such a tree structure. 



Proposition 6.1. 

In the case where all wave numbers equal kq = k\ = • • • = n n , we have (IIT) n = (TII) n = 0 
where n is the number of sub-domains. 

Proof: 

First of all note that (TII) n = n(IIT) ra n, so we only need to prove the result for (IIT) n . 
We start by simply writing down the explicit expression of the operator (IIT) n . For any 
u,C6 H(X), setting := R],(Uj) and 0 ^. := R/(,(t)j), we have 

n 

[nT( U ),t>] : =E E E 

*o=0 

*2 7^*0 
n 

i(nT) 2 (u),t)] :=e X X X 

*o—0 

* 27^*0 * 37^*1 

In the expression above we have i -2 / io and is i\ due to the very definition of the operator 
T. Applying recursively the formulas derived above for FIT finally leads to the following 
explicit expression for (IIT) n , 
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( 21 ) 


i(nT)“(u),»] ==£ v v £ E. E 

* 0=0 iiE^ig * 26 ^ *3€^»2 MG^tg Ls+lG^in 

* 27^*0 *3^*1 * 37^*2 *»+l^*n-l 

,i2Q J ^‘U,*3 ' ' ' ^*0,nlri 0l ii 

The sum in the expression above is taken over all the sequence of indices io, h, ■ ■ ■, i n +i 
satisfying the constraints € d^i k _ 1 (which implies in particular that / ik-i) and i*. +1 / 
ik~ 1 for all k = 1... ra. Each %k is the index of the sub-domain and Hj fc+1 is adjacent 
to Qi k since ik+ 1 € J^i k , hence those sequences io, h, ■ ■ ■, i n +i are actually chains of length 
exactly ra + 1 of the tree T. But since T only admits n elements and is a tree, it does not 
contain any such chain. This implies that the summation in (21) is taken over an empty set. 
Hence [(nT) n (u), oj = 0, and since u, D were chosen arbitrarily, this finally implies (nT) n = 0. 
□ 

The spectrum of the local multi-trace operator can now easily be deduced from what precedes, 
in the case where all wave numbers equal. Note that the next result states equality and not 
just inclusion. 

Theorem 6.1. 

Assume that all wave numbers are equal i.e. kq = k\ = 
point spectrum of A — all i.e. the set of its eigenvalues 
coincides with the point spectrum, and it is given by 

©(A - an) = 6 P (A - an) = {+ \/l + a 2 ,-yjl + a 2 }. 


• • • = K n . Let & P (A — an) refer to the 
. Then the spectrum of this operators 



Proof: 

The result is clear if a = 0 since (Id + A)/2 is a projector, so for the remaining of this 
proof, we will assume that a/0. Taking the square of the above operator, and using Lemma 
6.2, yields (A — an) 2 = (l + a 2 )Id — a(Tn + nT). Then it is a direct consequence of Corollary 
6.1 and Proposition 6.1 that the operator Tn + nT is nilpotent. Hence according to [21, 
Thm.10.13] and the spectral mapping theorem [21, Thm.10.28], we have 6(Tn + nT) = {0}. 
This also shows that ©((A — an) 2 ) = {1 + a 2 }, hence applying once again the spectral 
mapping theorem, we finally conclude that ©(A — an) C {+Vl + a 2 , — V 1 + a 2 }. 

Denote for a moment /(A) := A 2 . Then clearly © p (Tn + nT) = {0} since Tn + nT is 
nilpotent. Moreover, according to the ’’point spectrum counterpart” of the spectral mapping 
theorem [21, Thm.10.33], we have /(© p (A — an)) = © p (/(A — an)) = {1 + a 2 }. Hence we 
conclude also that © P (A — an) C {+VT+~c?, —a/ 1 + « 2 }- If we can prove that the previous 
inclusion is actually an equality, then the proof will be finished. It suffices to show that, if 
A € © p (A — an), then we also have —A € © P (A — an). 

Take an eigenvector u € H(E) \ {0} of A — an associated to the eigenvalue A. Since A 2 = 

n 2 = id, we have A(An - nA) = -(An - nA)A and n(An - nA) = -(An - nA)n. As 
a consequence (A — an)(An — nA)u = —(An — nA)(A — an)u = — A(An — nA)u. Hence, 
—A is an eigenvalue of A — an if A is, provided that (An — nA)u / 0. Observe that, since 
(A —an)u = Au and A 2 = n 2 = Id, we have nAu = au + Anu and Anu = (l/a)u— (A/a)Au. 
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Summing these two identities yields 


(IIA — AII)u = (a — l/a)u + — (A + all)u. 

a 

If we can prove that 0 is not an eigenvalue of (a — l/a)Id + ^(A + all), this will show 
that (IIA — AII)u V 0. From the first part of the proof, we know that the spectrum of 
A + all is included in {+\/l + a 2 , — Vl + a 2 }. Besides A equals +Vl + a 2 or — Vl + a 2 . As 
a consequence the spectrum of (a — l/a)Id + ^(A + all) only contains the values 


a-± —\/l + a 2 

a a 


a 2 — 1 ± (1 + a 2 ) 
a 


2a or — 2/a 


Since a € C \ {0}, we have 2a V 0 and 2/a VO. As a consequence the spectrum of 
(a — l/a)Id + ^(A + all) does not contain 0 so (IIA — AII)u V 0 necessarily. □ 


The previous theorem can be reformulated as follows. 

Corollary 6.2. 

Assume that all wave numbers are equal k,q = = ■ ■ ■ = K n . Then for any pair of complex 

numbers a, f5 € C the operator A + all + /3Id is invertible if and only if (3 2 — a 2 V 1. 

Theorem 6.1 also leads directly to explicit expression for the spectrum of the local multi-trace 
operator. Thanks to Fredholm theory, this implies a well-posedness result. 

Corollary 6.3. 

For any a £ C \ {0}, the operator L a := (1 — a)(A — Id) + a(A — II) is invertible. Moreover, 
in the case where all wave numbers are equal kq = = ■ ■ ■ = n n , its spectrum equals its point 

spectrum and ©(L Q ) = © p (L a ) = {—1 + a — Vl + a 2 , — 1 + a + Vl + a 2 }. 

Proof: 

Assume that the wave numbers Ko,Ki,...,n n are arbitrary elements of (0,+oo), and 
consider any a £ C \ {0}. Let A* refer to the operator defined in the same manner as A but 
with wave numbers all equal to k* = i. Then the operator A—A* : BI(S) —>• H(£) is compact as 
it only involves integral operators with regular kernels, see [22, Lemma 3.9.8]. Then Theorem 
6.1 shows that (1 — a)(A*— Id) + a(A* — II) is invertible, since its eigenvalues —l + a± Vl + a 2 
differ from 0 as a V 0. Hence (1 — a) (A — Id) + a(A — n) is a compact perturbation of an 
isomorphism. According to Fredholm-Riesz-Schauder theory (see [16, Chap.2] for example), 
this operator is invertible if and only if it is one-to-one. Since it is injective according to 
Proposition 4.1, we finally conclude that (1 — a)(A — Id) + a(A — n) is an isomorphism. The 
second statement above concerning the spectrum is a trivial consequence of Theorem 6.1. □ 


In the case where wave numbers take arbitrary values the spectrum is not reduced to — l + a± 
Vl + a 2 anymore. However a difference of wave numbers only induces compact perturbation 
of integral operators so that, in the general case, this result still indicates the location of 
accumulation points of the spectrum. 

Corollary 6.4. 

For any a £ C\{0}, set L a := (1 — a)(A—Id)+a(A—n). Then any element A £ 0(Lq,)\{— 1 + 
a + Vl + a 2 , —1 + a — Vl + a 2 } is an isolated eigenvalue with dim( ker(L a — Aid)) < +oo. 
Moreover the two values —l + a± Vl + a 2 are the only possible accumulation points of &( L a ). 
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Proof: 

Denote := — 1 + a± \/l + a 2 , and set £(A) := L a — Aid. Then A i->- £(A) is an analytic 
operator pencil, and it is Fredholm valued for A 7 - Indeed take any A £ C \ {/r+, p,~}, 

and define the operator in the same manner as L a except that all wave numbers are taken 
equal to kq- The operator £'(A) := L' a — Aid is invertible according to Corollary 6.3, and 
£'(A) — £(A) is compact. As a consequence, £(A) is a compact perturbation of an invertible 
operator, so it is Fredholm with index 0 and admits finite dimensional kernel. 

Further, since L a is a bounded operator, £(A) is invertible for A > ||L a ||, where || || refers 
here to the norm naturally associated to continuous operators mapping H(E) into itself. As a 
consequence, we can apply Fredholm analytic theorem (see Appendix A in [12]) which shows 
that £(A) is invertible in C\ {/r+, p~} except for a countable set of isolated values. Moreover, 
we have just seen that all these values can only lie in the disc of center 0 and radius ||L Q ||. □ 

7 Numerical evidences 

In this section, we present a series of numerical results confirming the conclusions presented 
previously. We consider 2-D scattering problems of the form (2) involving three domains. As 
regards discretization, we consider a uniform paneling T, h ~ E which induces a mesh for each 
of the sub-domains Tj ~ Tj, Tj C E ft . The discrete spaces H/ l (E) are constructed on these 
meshes by means of Pi-Lagrange shape functions for both Dirichlet and Neumann traces 

%(E) = { (u h j,p h j ) j= 0 , 1,2 such that 

( 22 ) 

Vj = 0,1,2,for all panel e C T^, Uj\ e ,Pj\ e € Pi(e) }. 

Denote the matrix associated to the Galerkin discretization of the local multi-trace for¬ 
mulation ( 12 ) by means of the discrete space ( 22 ), and let us denote the matrix obtained 
by Galerkin discretization of the bilinear form (u, t») 1 —> [u, t>|. We shall focus our attention 
on the spectrum of the matrix (M/,) _ 1 B/j that may be considered as an approximation of the 
continuous operator associated to Formulation (12). 


Do 


Figure 1: First geometrical configuration 

All computations have been achieved on a laptop equipped with a 2-core Intel i7-3520M 
processor at 2.9GHz with 4 GB of RAM. Meshes have been generated using Gmsh [9] (see 
also the website http://geuz.org/gmsh/). For computation of eigenvalues we relied on the 
Arpack+-1- library (see http: //www. ime .unicamp .br/~chico/arpack++/). 
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Figure 1 represents a first geometrical configuration. The boundary of fio is a unit square 
centered at 0, and the boundary of is a circle of radius 0.5 centered at 0. Figure 2 represents 
the spectrum of for a mesh width h = 0.05 and Kq = K\ = = 1- 




Figure 2: Two clusters of eigenvalues in the case of no contrast of wave numbers (left) and a 
zoom on the cluster of positive eigenvalues (right). We took a = 1 and kq = n\ = K 2 = 1. 


The spectrum clearly takes the form of two clusters centered at the values zb\/2 = — 1 + 
a ± n/1 + a 2 for a = 1. Figure 3 shows the spectrum of the same matrix, with the same 
geometrical configuration, but with a = 0.5 and a = —0.25. The formula —1 + a zb \/l + a 2 
yields the values 0.61803 and —1.6180 for a = 0.5 (up to 5 digits), and —0.21922 and —2.2808 
for a = —0.25, which is consistent with our theory. 




Figure 3: Spectrum of the local multi-trace operator with ko = Hi = = 1 and a = 0.5 

(left) or a = —0.25 (right) 


Next, in Figure 4, we consider the case a = 1, but wave numbers differ taking the values 
kq = 1, Ki = 5 and K 2 = 2. The mesh width remains h = 0.05. Although the eigenvalues are 
not clustered anymore, they are more densely grouped around zb\/2 suggesting that these are 
the only two accumulation points of the spectrum of the continuous operator. 
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Figure 4: Spectrum of the local multi-trace operator in the case kq = 1, = 5 and K 2 = 2 

with a = 1. 

For the next series of figures, we consider the same scattering problem, but in a different 
geometrical configuration. The new configuration is depicted in the picture below: there 
are two square scatterer separated by a thin gap of width delta. In the theory we have 
presented, we needed to assume that there is no junction point i.e. points where at least 
three sub-domains abut. We wish to test the robustness, with respect to this assumption, of 
the theoretical formulas obtained. 


( —1, +0.5) 


(1, +0.5) 






In Figure 5, we consider kq = «i = k 2 and a = 1. For a fixed strictly positive value of 5 > 0, 
the eigenvalues are clustered around ±\/2- Each of the figures below represent a zoom at the 
cluster centered at +y/2 for various values of 5. 




















00 




( a ) ( b ) (c) 

Figure 5: Spectrum of the local multi-trace operator for kq = n\ = K 2 = 1, a = 1 and three 
different values of gap: S = 0.1 (left) 5 = 0.01 (center) and 5 = 0.001 (right) 
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The cluster is more and more scattered as the gap closes, suggesting that the assumption that 
there is no junction point is mandatory, in spite of quadrature rules being less reliable as 6 is 
close to 0. In the last picture below, we examine the case where the gap is closed i.e. (5 = 0, 
which corresponds to the presence of a junction point in the geometry. 




fl 0 


Figure 6: Spectrum of the local multi-trace operator (left) for kq = k\ = K 2 = 1, a = 1 in 
the presence of junction points in the geometry (right). 


Unfortunately, geometrical configurations involving junction points are not covered by the 
theory of the present article. However the result of Figure 6 suggests that, although the 
eigenvalues are not anymore closely clustered around the values — 1 + a ± \J\ + a 2 , these 
two theoretical points remain accumulation points of the spectrum of the local multi-trace 
operator. 

Acknowledgement The author would like to thank V.Dolean and M.J.Gander for fruitful 
discussions. 
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